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SUMMARY 


The endochronic theory of plasticity originated by Valanis has been 
applied to study the axially symmetric motion of circular cylindrical thick 
shells subjected to an arbitrary pressure transient applied at its inner 
surface. The constitutive equations for the thick shells have been obtained. 
The governing equations are then solved by means of the nearcharacteristics 
method . 


INTRODUCTION 


The problem of dynamic plastic response of shells has received consider- 
able attention in recent years. Most of the published works are based on the 
flow theory of plasticity and usually limited to Isotropic linear work- 
hardening materials. Theoretically, the flow theory is based oh the existence 
of an initial yield surface coupled with an assumed hardening rule to obtain 
subsequent yield surfaces; an extensive bookkeeping is necessary to trace the 
evolution of the yield surface which changes as deformation progresses. The, 
analysis of inelastic responses of the bodies is therefore complicated by path 
dependence and the yield condition, which introduces different governing 
equations in the distinct regions - elastic and Inelastic. Valanis (ref. 1) 
presented a new theory of plasticity termed endochronic theory, which com- 
pletely abandoned the concept of a yield surface and its subsequent hartfening 
rule. 

The endochronic theory of plasticity is based on thermod 3 mamic theory of 
internal variables and conforms to experimentally observed material behavior. 
The basis of the endochronic theory is the assumption that the current state 
of stress is a functional of the entire history of deformation. The influence 
of past deformation on the current stress is measured in terms of a mono- 
tonically increasing time scale of strain-defined (ref, 1) or stress-defined 
(ref. 2) endochronic time. This theory has been applied to give analytic pre- 
dictions for the quasi-static mechanical response of engineering materials 
(metallic (ref. 3) and non-metallic (ref. 4)), the dynamic response of a 
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thin-walled tube subjected to a combined longitudinal and torsional step 
loading (refs. 5,6), and the dynamic plastic response of circular cylindrical 
thin shells (refs. 7, 8). It has been shown that the theory does indeed have 
the capability of explaining the observed phenomena quantitatively with 
sufficient accuracy. 

In this paper, the endochronic theory is applied to thick axially- 
S 3 rmmetric cylindrical shell subjected to dynamic loading. The governing equa- 
tions are then solved by the method of nearcharacteristics. 


FORMULATION OF THE PROBLEM 


Consider a circular cylindrical thick shell with mean radius R and 
thickness H. For the axisymmetric motion of shell, the stress and strain 
states are 
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where ^ is the Cauchy stress tensor, s is small strain tensor, s and ^ are 
the deviatoric stress and strain tensors, respectively, and subscripts x, r, 

0 refer to the components in longitudinal, radial, and circumferential 
directions, respectively, Let U and W denote the displacements in the axial 
and radial directions respectively at time t of the cross section a distance 
X from a referenee section, and u and w are the correspohding velocity 
components. The equation of motion in the x and r directions have the 
following form; 
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where p is the density . 
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The strain-displacement relations and the Gorresponding compatibility con- 
ditions are 
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For isotropic material under isothermal condition with the assumption of 
elastic hydrostatic response, the constitutive equations in the endochronic 
theory can be found from reference 1 as follows: 
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where a, 3, K^, K 2 are the material parameters, p is shear modulus, K is bulk 
modulus, kk, iZ, and ij are subscripts denoting coordinates, d^ is the 
endochronic time measure with the restriction that + K 2/3 > 0, K 2 t. 0, 
and Ki and K£ may not both be zero. From the definition of s and ^ considered 
in this problem, it is possible to express the time measure approximately as 


gd^ = + 





. 1/2 


"^9 


( 12 ) 


where Bi = E^/Oq, is the asymptotic slope of the uniaxial stress-strain 
curve for large strain, Oq is the intercept of this slope with the stress 
axis, and the positive sign holds for straining while the negative sign is for 
unstraining of de 0 . Using (12i and equations (9), (10) and the compatibility 
conditions (5) to ( 8 ) results in the following: 


3a 


X 


3t 






3a 3a 3a„ 

X , r , 9 _ 

3t " 3t 3t “ ^2 


(13) 

(14) 


451 



9a 9a 9a. » . 

_x . __r . e _ „ 9w ^ 

9t ^ 9t + 3t 9x 9r ^3 


9a „ 

xr _ ^ 

9t ^ 9r 


9w ^ 
^ 9^ = ^ 


(15) 

(16) 


where 


^2 = 


- a. (2a -a -a.) 

- i 1 X r 9 w 

+ 2 1 + e? r 

■„ - 1 

^ 2 1 + 3 C 


w 

r 


aj . 3K f 

- 2 xr w 

^4 + 1 + e? r 


“l = 




de 


Qj 


^de 
1 

,de. 


1/2 




- 6j_a 


and E is elastic modulus, v is Poisson's ratio. Equations (3), (4), and (13) 
to (16) are the fundamental equations of the problem considered here. 


NEARCHARACTERISTIC SOLUTION 


The governing equations presented in the previous section together with 
the auxiliary equations can now be written in matrix form as follows: 
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and 


{X} =■ 
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The above set of equations is of hyperbolic type j the conventional 
bicharacteristic method would be very tedious for six dependent variables . 

Using the method of nearcharacteristics first proposed by Sauer (refi 9), we look 
for characteristlc-llke lines in the coordinate planes along which the solu- 
tion can be extended. (Sauer called these lines nearcharacteristics.) The 
formulation and numerical technique in the nearcharacteristics resembles the 
one-dimensional approach except that those partial derivatives which do not 
lie in the plane of interest are considered of zeroth order in that particular, 
calculation. For example, when the bicharacteristics in the x-t plane are of 
interest, then those terms in [A] containing partial derivative in r-direction 
are combined with terms in {B} in equation (17) . Now following the same 
procedures as described in reference 8 for one-diniensional case, the near- 
characteristics in the x-t and r-t planes, respectively, are obtained as follows 
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The near character is tics obtained here indicate that there are two character- 
istic cones existing in the present analysis; one of them (eq. (19)) corre- 
sponds to the longitudinal wave propagation while the other (eq. (20)) 
corresponds to shear wave. They are right circular cones with their center 
lines perpendicular to the x-r plane as shown in figure 1. This is an expected 
result, because the governing equations have constant coefficients for the 
highest order terms. There are no convected terms appearing in the present 
analysis. The compatibility equations along the nearcharacteristics can be 
found in the same way as in the one-dimensional case. In the x-t plane, we 
have; 
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Note that each set of the above characteristics lies entirely in planes paral- 
lel to one of coordinate planes. Equations (21) to (28) have the appearance of 
a one-dimensional method of characteristics formulation except that they con- 
tain the partial derivative terms in the other coordinate direction. The 
nearcharacteristics equation derived here can be solved nximerically by the one- 
dimensional teehnlque. Two independent solutions can be obtained, each corre- 
sponding to one of the coordinate planes. 

NUMERICAL EXAMPLE 


Consider a central segment of the Clinch River Breeder Reactor steam gen- 
erator flow shroud with length 21 = 1.0668 m, mean radius 0.47 m, and thickness 
0.0127 m. The material is 2.25 Cr-1 Mo at 756 K. The pressure input function 
was generated by the hydrodynamics module (ref. 10). A constant volume, step 
pressure pulse of 13.79 MPa was taken as the source pressure p at the center. 
This is typical of the maxima observed in large sodium-water reaction experi- 
ments during the transient period. Since the pressure loading was supposed to 
be S3amnetrie with respect to the mid-span, only half-length of the shell needed 
to be considered here. The boundary conditions for the example are shown in 
figure 2- as follows: 
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It has been shown in reference 11 that the two independent solutions, each 
based on one coordinate plane, are numerically unstable liidiile a calculation 
method obtained by averaging the above mentioned independent solution yields a 
stable solution. In view of equations (21) , (22) and the boundary conditions 
(29), it appears that the neareharacteristics equations in x-t plane are not a 
proper choice at r = 0 and r = H because Oj- are being prescribed there. There- 
fore a combination technique is proposed here: on the boundaries r = 0 and 
r = H the solutions are obtained from r-t plane nearcharacteristics equations 
while at other points the solutions are obtained from the x-t plane. The 
ntimerical results here show that this leads to a stable solution. The advan- 
tage of this technique over the averaging method is a tremendous saving in 
computation time. The resulting pressure history at the midspan (x = 0) of the 
middle surface of the shell is shown in figure 3. The resultant dynamic 
response of radial displacement and velocity as a function of time for the same 
center point of the shell is also shown in the figure. In figure 4, shell 
displacement profiles are shown for several times. 

CONCLUDING REMARKS 


The endochronlc theory of plasticity originated by Valanis has been ap- 
plied to study the axially symmetric motion of circular cylindrical thick 
shells subjected to an arbitrary pressure transient applied at its inner 
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surface. The constitutive equations for the thick shells have been obtained. 
The governing equations are then solved by means of the nearcharacteristics 
method. It has been shown that a stable solution can be obtained by treating 
the radial boundaries in one coordinate plane while at other points the solu- 
tions obtain from the other coordinate plane, 

REFERENCES 


1. Valanis, K. C. : A Theory of Viscoplasticity without a Yield Surface, Part 

I, General Theory. Arch. Mech. Stosow. , vol. 23, 1971, pp. 517-534. 

2. Valanis, K, C. ; On the Foundations of the Endochronic Theory of Visco- 

plasticity* Arch. Mech* Stosow., vol. 27, 1975, pp. 857-868. 

3. Valanis, K. C. ; A Theory of Viscoplasticity without a Yield Surface, Part 

II, Application to the Mechanical Behavior of Metals. Arch. Mech, 
Stosow., vol. 23, 1971, pp. 535-551. 

4. Bazant, Z. P. ; A New Approach to Inelasticity and Failure of Concrete, 

Sand and Rock-Endochronic Theory. Proceeding of 11th Annual Meeting, 

Soc. Eng. Sci. , Duke Univ. , 1974. 

5. Wu, H. C. ; and Lin, H. C, : Plastic Waves in a Thin-Walled Tube under 

Combined Longitudinal and Torsional Loads. 10th Anniversary Meeting, 

Soc. Eng. Sci., North Carolina State Univ., Nov. 1973. 

6. Wu, H. C. ; and Lin, H. C. : Combined Plastic Waves in a Thin-Walled Tube* 

Int* J. Solids Struct., vol. 10, 1974, pp. 903-917. 

7. Lin, H. C. : Dynamic Plastic Deformation of Shell - An Endochronic Solution. 

Trans. Am. Nucl* Soc., vol* 22, 1975, pp. 567-568. 

8. Lin, H. C. ; Dynamic Plastic Deformation of Axi-Symmetric Circular 

Cylindrical Shells. Nucl. Eng. Design, vol. 35, 1975, pp. 283-293* 

9. Sauer, R. : Differenzenverfahren fur hyperbolishe Anfangswertprobleme bei 

mehr als zwei unabhangigen Veranderlichen mit Hllfe von Nebencharak- 
teristiken. Numerische Mathematik, vol, 5, 1963, pp. 55-67. 

10. Shin, Y. W. ; and Valentin, R. A.: Dynamic Structural Loads Produced by 

Large Sodium-Water Reaction in LMFBR Steam Generators. 3rd SMiRT, 

London, UK, 1975, Paper E4/5*. 

11* Shin, Y. W. ; and Kot, C. A.: Two-Dimensional Fluid-Hammer Analysis by the 

Method of Nearcharacteristies. Argonne National Laboratory, ANL-75-21, 
1975. 


456 




PRESSURE 

4 

0 

X = 0. 


m 

0.035 
0.030 
0.025 

RADIAL 0.020 
DISPLACEMENT 

0.015 

0.010 

0.005 

0.000 

0 0.1 0.2 0.3 0.4 0.5 0.6 

X. m 



DISPLACEMENT 



RADIAL 

VELOCITY 


1000 

TIME, |is 


Figure 3.- Radial displacement velocity, pressure history at 


Figure 4.- Radial displacement profiles. 



